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SUMMARY 


Various investigations of single -degree-of -freedom pitching oscil- 
lations of a wing in potential flow have heen made. However, corre- 
spondlng studies of single-degree-of -freedom flutter of a control 
surface have not "been made. The present paper demonstrates "by theo- 
retical calculations that single-degree-of-freedom control-surface 
flutter is possible. The effects of structural damping, aerodynamic 
balance, eixis of rotation, and compressibility eire included. 

The mechanism of instability presented here is based on potential- 
flow theory and the results of such a study are not directly applicable 
to separated flowj in addition, certain practical limitations are 
discussed. 


INTRODUCTION 


Flutter of a control siirface alone, that is, an unstable oscilla- 
tion in which the control surface is effectively rigid and nondeformable 
and rotates around its hinge line, has been encountered on aircraft in 
flight. This phenomenon, as well as stall flutter, is thought to be 
associated with separated or nonpotential flow, although the explanation 
is not entirely clear. Flutter of a control surface in one degree of 
freedom is possible in potentlaL flow, at least theoretically, and forms' 
the siibject of the present paper. The development of this subject may 
possibly contribute to the understanding of the nonpotential-flow 
problem . 

Theoretical investigations of several types of undamped single- 
degree-of-freedom oscillations in potential flow have been reported 
recently. For instance, Smilg (reference l) has reported a study of 
the tindamped oscillations of a pitching wing in incompressible flow and 
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this study has iDeen extended ty Runyan in reference 2 to include the 
effect of compressibility. Cvmningham (reference 3) has studied the 
undamped bending oscillations of a swept wing in subsonic flow. The 
undamped pitching oscillations of a wing in supersonic flow have been 
studied hy a numher of investigators, for Instance, Garrick and Ruhinow 
(reference h). Temple and Jahn (reference 5)^ and Watkins (reference 6). 

The aforementioned papers concern single -degree-of -freedom pitching 
of a wing or the bending of a swept wing. Some examples of control- 
surface flutter in one degree of freedom are given in reference 7 along 
with a discussion of some possible causes. Apparently no calculation 
studies exist demonstrating the conditions for the existence of single- 
degree-of-freedom oscillations of a control surface. Accordingly, the 
purposes of this piaper are (l) to show the conditions for the existence 
of slngle-degree-of -freedom oscillations of a control surface based on 
potential-flow theory and (2) to examine and present the effect of various 
parameters such as Mach number, fluid density (or altitude), location of 
axis of rotation, structural damping, and aerodynamic balance on this type 
of oscillatory instability. 

Certain limitations which may affect practical consideration of the 
results presented in this paper should be noted: (l) the calculations 

are based on two-dimensional aerodynamic coefficients, end the effect 
of aspect ratio, which could be appreciable, is not taken into account, 

(2) the influence of additional degrees of freedom, which will always 
be present in actual configurations,- is not considered, and (3) the 
effect of finite thickness and shape of the control surface is not 
considered, since the control surface has been replaced in theory by 
its mean camber line. In addition to these limitations, the results 
of this study are not directly applicable to separated flow, since the 
analysis was based on potential -flow theory. Nevertheless, this approach 
should provide a basis for coD5)arison and study and a guide for further 
investigation. 


SYMBOIS 


wing half -chord, feet 

location of hi n ge axis of control surface without aerodynamic 
balance with respect to midchord point, based on half -chord 
and positive rearward 

location of leading edge of control stirface with aerodynamic 
balance with respect to midchord point, based on half -chord 
and positive rearward 
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Cp torsional stiffness of control surface about hinge line c 

Cj^ function of ratio of aileron chord to wing chord 


e location of hinge axis of control surface with aerodynamic 

balance with respect to midchord point, based on half -chord 
and positive rearward 


F and G 
gp 

^bp 

"P 


transcendental functions of k for oscillating plane flow 
structural damping coefficient ~ Logarithmic decrement) 

out-of -phase aerodynamic moment coefficient 

total out -of -phase moment coefficient 

mass moment of Inertia of control surface about hinge line 
per unit length 

^RR aerodynamic moment coefficients for a control surface In 
compressible flow 


k 


reduced-frequency pa3ameter (bio/v) 


2 nondimensional. distance from control -surface axis of rotation 

to leading edge of control surface (e - c') 

M Mach niunber 


M. 


total complex aerodynamic moment on control surface 


^p 

%p' 

T 

n 

V 

P 


aerodynamic coefficients for a control surface in supersonic 
flow 

total inphase moment coefficient 

part of Inptiase aerodynamic moment on control surface 

control -surface coefficients defined in references 8 and 9 
flutter velocity, feet per second 

control- surface rotation, measured from wing chord line 
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fluid density, slugs per culic foot 
flutter frequency, radians per second 

natirral frequency of control surface about hinge eixis, radians 
per second 


AMiySIS 


In this section, the ejqjression for the equilibrium of moments on 
a control surface oscillating about its hinge line is given, followed 
by a discussion of the aerodynamic moment for (l) a control surface 
oscillating in Incompressible flow, (2) a control surface, having aero- 
dynamic balance, oscillating in incompressible flow, and (3) a control 

surface oscillating in compressible flow. 

% 

The equation representing the equilibrium of moments per unit 
length on a control surface oscillating about the hinge line c is 


Ipp + (1 + igp)Cp3 = Mp 


( 1 ) 


where Ip is the mass moment of Inertia per unit length of the control 

svirface aboxrt its hinge line c, P is the angular deflection of the 
control surface measured from the chord line, gp is a structural 
damping coefficient, and Cp is the torsional stiffness of the control 
system about c. The quantity Mp is a complex aerodynamic moment, 
which is a function, in part, of the displacement 3, velocity 3, 
acceleration 3^ reduced-frequency parameter k, Mach number M, and 
location of hinge eixis c. T^ complex notation has been in use in 
flutter work for many yeeirs because of its convenience, since it contains 
both the magnitude of the forces or moments and the phase with respect 
to the displeicement . The escpression for the aerodynamic moment Mp 
is given in reference 8 for a control surface which is hinged at its 
leading edge and oscillating in incompressible two-dimensional flow. 

An e3q>ression for the aerodynamic moment Mp for an aerodynamlcally 

balanced control surface is given in reference 9* Closed esspressions 
for Mp for compressible fluid do not exist and values of the aero- 
dynamic moment must be obtained by use of tables as given, for 

exan 5 )le, in references 10 and 11 for subsonic flow (M < l) and in 
reference 12 for supersonic flow (M > l) . 
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Since equation (l) is complex, it can "be separated into two com- ' 
ponents as follows: 

° (2) 

"where Rj^p represents the inphase (real) moment coefficient and 
represents the out-of -phase moment coefficient for the_control sxirface. 
Equation (2) in5)lies that both components, I-^p and must vanish 

simultaneously for equilihrium. 

The vanishing of the out-of -phase moment corresponds to a 

borderline condition between damped and undamped oscillations as was 
more fully discussed in reference 2 . The flutter frequency may then 
be determined from the inphase (real) moment equation. This Inphase 
moment equation = 0 may be put in more convenient form as follows: 




( 3 ) 


where R^^^ ’ 


is proportional to the Inphase aerodynamic moment on the 


control surface. In the sxibsequent sections, expressions for 
and R^p' will be given for the various cases. 



Equations for Incompressible Flow 

Control siirface wlthoiib eterodynamlc balance .- The expression for 
the aerodynamic moment on a harmonically oscillating control surface 
in Incompressible flow having no aerodynamic balance may be derived from 
reference 8 as 


Mp = -pV 


i v2(t5 - TljTio)p - ^ vbWliTu - i 




(W 
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■where the T coefficients axe functions only of the hinge position c 
and are defined in reference 9* Substitution of equation (it-) into 
equation (l) aind separation into the two conrponents lead to eqiiatlon (2) 
where 


®hp = 


'b. + 

„2 


T^Tio) - 


*^^12/^11 ^ ^ ^ 

2n \2rt k 


► k 


:)+ A 

- 1 

7 JTpb^ 



= 0 


(5) 


and 




T12P1O 2G 


2 it 


2it 


^ ^ 2F - 


2jr 


2 ^it-^n 


gp 


npb 


= 0 


( 6 ) 


The Inphase (real) moment may be expressed in the form given by 

equation (3) where 


Rbp’ = 


^12 / ^11 ^ 

2ir \2n k 



(7) 


A numerical or a graphical method is necessary for the solution of 
the out-of-phase moment I-pp = 0 since the functions F and G are 

transcendental functions of l/k. Once the combination of l/k and 
F and G that satisfies equation (6) has been determined, the frequency 
and fl'utter speed can be obtained from equation (3) in which the value 
of Epp' is taken for the same valtie of l/k that satisfied equation (6). 

Control surface with aerodynamic balance .- The force and moment 
on a harmonically oscillating control surface having aerodynamic 
balance ha've been derived in reference 9 tor incon^iressible flow. A 
control sttrface having aerodynamic balance is represented in refer- 
ence 9 as follows: 
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where the dimensions c, c% e, and 1 are referred to the half- 
chord h . 

The out-of -phase (imaginary) moment coefficient may he derived 
from reference 9 as 



= 0 ( 8 ) 

and the freguency is determined from the inphase (real) moment expressed 
as equation ( 3 ) where 

^3' ' - ^(-..''3 - *%) + ^(’’18 + + *^28) - 

4(Ti 2 - 21l2o)[^(Tii - - i(Tio - JTsi)^] (9) 

The T coefficients are given in the appendix of reference 9 and 
all T's are functions of c. One exception that needs separate con- 
sideration is the term Tg^ which is dlsciissed in reference 9* For the 

case (over idealized) of a sharp vertical step, the term T 20 becomes 

infinite; however, for practical configurations certain approximations 
can be made which permit a finite-value determination of this term. 

One such simple approximation is given in reference (9) as 

c’ = ^e + 3c) 


Equations for Compressible Flow 

Closed expressions for aerodynamic moment and lift on oscillating 
wings and control surfaces in compressible flow do not exist. Aero- 
dynamic moment and lift coefficients have been calculated by severail 
investigators for various ranges of Mach number and reduced frequency 
and are us\ially presented in the form of tables. In reference 10, 
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tables of aerodynamic coefficients are presented for one aileron-to- 
ving-chord ratio of 15 percent (c = O. 7 ) for M = 0.7- These tables 
are extended in reference 11 to include several aileron-to-wing-chord 
ratios of 2k, 33^ and k2 percent. 

A notation differing from that used in this paper is used in 
references 10 and 11 -where the inphase (real) moment coefficient of 
a control surface is denoted hy and the damping (out-of -phase) 

moment is denoted hy The relation between the coefficients of 

this paper and those of references 10 and 11 is given as 

= 0 ( 10 ) 

U t _ 

where Cj^ is a function of the ratio of the aileron chord to the wing 
chord and is given in reference 11, page I 6 . 

In order to obtain a solution, it is necessary to plot against 

k and determine the -value of k at which = 0. The frequency is 

then gl-ven for the same valvie of l/k by 



The resxilts of calc-ulatlons of aerodynamic coefficients for a 
control surface oscillating in supersonic flow ha-ve been reported in 
reference 12. The aerodynamic moment on the control 'surface is gi-ven 
as 


Mp 


-lipb^v^k^P^N^ + IRgj 


( 12 ) 


where is the Irphase moment coefficient and Ng is the out-of- 

phase moment coefficient. In order to obtain a solution for the out- 
of-phase moment, Rg = 0, the value of l/k at which Rg = 0 must be 
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determined. The frequency is then given for the same value of l/k 

■by 



1 - 


5 I, 




( 13 ) 


DISCUSSION OF EESULTS 


In this section, the results of calculations based on the expres- 
sions previously presented are given. The oscillation of a control 
surface in Inconipressihle flow with zero structural damping is discussed 
and the effect of structural damping is Introduced. A discussion of 
the effect of aerodynamic balance and the effect of compressibility 
follows . 

Flutter of control surface with' zero structural damping in incom- 

piesslble flow .- Calculations based on equations ( 5 ) and (6) where 
gp = 0 have been made for various locations of the control-surface 
axis of rotation. The results of the computations are presented in 
figures 1 and 2 for ratios of the control- surface chord to wing chord 
of 10 percent, 30 percent, 50 percent, and 100 percent (c = 0.8, 0.i^, 

0, and -l) . 


In figiire 1 the flutter- speed parameter v/bo:p is plotted against 

an inertia parameter Ip/jtpb^. The curves represent the transition, 
for a given Mach number, from a damped or stable condition to an undamped 
or unstable condition. The stable region is below or to the left of 
a ciirve, and the unstable region is above or to the rl^t of a curve. ■ 
For small values of the inertia parameter the control surface would be 
stable and, if the inertia parameter is increased, for instance by an 
increase in altitude, a val\ie equal to the vertical asyn^rtote could be 

Ir 

reached or exceeded. In figure l(a) (c = 0.8) this value is — = 7*56 

Jtpb^ 

which increases to 577*7 (fig* l(*i)) for c = -1.0, that is, for a 
wing oscillatii^ about its leading edge. The flutter speed for values 

of Ipy^rtph^ equal to the asymptote is infinite, but a slight increase 

in Ipy/npb^ would result in a very rapid decrease in the flutter speed. 

For values of Ip/itpb^ approaching infinity, v/bci^ approaches 

a value equal to v/bu) at which the oscillation occurs. It should be 
noted that the oscillation occurs at a constant value of v/bco, which 
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Tiip nnFi -btnatj no matter wiiat the value of the inertia parameter^ the wave 
len^h of the oscillating wake will he the same in the border condition. 

The frequency ratio (cu/ci:p)^ is plotted against the inertia' param- 
eter Ip/jtph^ in figure 2 for the same axls-of -rotation locations as 
for figure 1. The unstable region is above or to the right of a curve. 
The vertical as y mptotes are the same as for figure 1, but the horizontal 
a sympt ote is unity. This fact indicates that^ for very large values of 
Ip^jtpb^, the frequency of oscillation approaches the natural (still air) 
frequency of the control siorface. 

In figure 3, the value of the reduced velocity l/k is plotted 
against the location of the axis of rotation of the control surface. 

Note that' the variation of l/k with c is not very large 

(25 < l/k < 35), but the values are considerably hi^er than experienced 

in flrtter anH approach values found in stability work. 


In figure 4, the minimum value of the inertia parameter (vertical 
asymptote) at which the oscillation could occur is plotted against 
location of control-surface axis of rotation. The value of the inertia 
parameter increases as the aileron-to-wing-chord ratio incieases. This 
plot is of particular significance with regard to the zero structural 
restraint case (cop = 0). If equation ( 3 ) is inverted to obtain 



-Rbp 



and ocp = 0 , there is obtained the relation 


Rbp' 


Jtpb^ 


m 


which Indicates that, if Ip/irpb^ is equal to ^ oscillation 

is possible. The frequency of oscillation is then a direct function of 
the velocity, and flutter can theoretically occur above zero airspeed 
as shown by the relation 


V 


bo) 

k 


where l/k is the value that satisfied Ij^p = 0. If the valtie of 
is exceeded by Ip/«pb\ the control surface is being operated past 
the transition point and in the xmstable region. 
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Flutter of control surface Td,th. structural damplrij:; In Incompressible 
flow.- As is often donej the type of danping force assumed is one that 
is Inphase with the velocity ^ hut proportional to the angular dis- 
placement P . 

!Fhe results of computations for one location of the axis of rota- 
tion c = 0.7 are presented in figures 5 and 6. In figure 5, the 
flutter-speed parameter v/hoip is plotted against the inertia parameter. 

As noted^ there is a very large effect of the damping parameter arid it 
may he that structural danping could he used for eliminating this type 
of oscillation. In figure 6, the frequency ratio (o)/aip)2 is plotted 
against the inertia parameter. It should he noted that flutter does 
not occur (contrary to the case for zero damping) at a constant value 
of the reduced speed v/hoi as the inertia parameter is varied. The 

values of v/hoi and Ip/irph^ for several values of the structural 
damping coefficient for c = 0.7 are as follows: 


gp 

l/k for - 

L_ 

- 20 

rtpb 

= 30 

itpb^ 

Ip 

mpb^ 

= 75 

mpb 

= 100 

' jrpb 

0 

31^.25 

34.25 

34.25 

34.25 

34.25 

.01 

38.0 

44.77 

56.1 

67.58 

7?. 60 

.02 

38.80 

48.36 

58.8 

71.50 

82.35 


Flutter of a control surface having aerodynamic balance and zero 
structural damping.- Calculations have been performed for two control 
surfaces by use of equations (8) and (9) for two locations of axis of 
rotation e hut one value of c, and the results are plotted in 
figures 7 and 8. In figure Z the flutter-speed coefficient is plotted 
against Ijj/jtph^ for c = O.i)- and e = 0.55 and e = O.75. The 
curves are similar to those obtained without aerodynamiic balance, except 
that the limiting value (vertical asymptote) appears to he at a larger 

value of the abscissa. In figure 8 the frequency ratio (oo/cup)^ is 
plotted against the inertia parameter Ip/nph^ for the same cases. 

The important facts to-be noted are (l) that aerodynamic balance 
did not eliminate single-degree-of- freedom instability of a control 
surface and (2) that the greater the amount of aerodynamic balance, the 
higher the limiting value of the inertia parameter at which the oscilla- 
tion could begin for this particular case. 
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FlTitrber of corrbrol surface in compreseible flov .- The results of 
the calculations based on the equations and tables discussed earlier 
are presented in figures 10^ and 11. In figure 9 the flutter-speed 
parameter is plotted agal.nrt the inertia parameter for M = 0, 0.7^ 
and 10 / 9 . The curves for M = O .7 and M = ^ are similar in form 

to the curve for M = 0, except that the stable range of the inertia 
parameter is greatly reduced as the Mach number is Increased. Another 
effect of increasing Mach nunber is to reduce the value of the reduced- 
velocity coefficient l/k^ aod this reduction is believed to be signifi- 
cant. In figure 10 the frequency ratio is plotted against the inertia 
parameter for three Mach nunber s M = 0, 0.7^ and I 0 / 9 . 

A significant plot is made when the asymptotic value of the inertia 
parameter is plotted against Mach number as shown in figure 11. The 
relation between the inertia parameter and Mach number appears to be 
linear^ and the region to the right and above the curve is the unstable 
range. This plot, as discussed earlier, is of partlciilar significance 
with regard to the control surface without elastic restraint (cup = O) 
and shows that an aileron that is stable in the low Mach number range 
could become unstable in the b1gb srobsonic or lower supersonic range. 

The asymptotic value of the inertia parameter at M = 0.7 for 
several positions of the control- surface hinge axis is included in 
figure 4. The effect of Mach number is again apparent in that increetsing 
Mach number considerably reduced the valtie of the inertia parameter for a 
given hinge position at which oscillation could occur. 


VARIOUS PRACTICAL LBDTATIOWS TO THEOKT 


In general, theoretical control- surface derivatives (either steady 
or unsteady) have not always been in good agreement with experiment. 

This difference is even more pronounced at the hi^er aircraft flight 
speeds and is partially due to breakdown and separation of the flow 
over the rear part of the wing. Since the aerodynamic coefficients 
have been derived on the basis of nonviscous, linearized, potential flow 
for simplified models, and the actual flow is viscous, nonlinear, and 
nonpotential, the differences between theory and experiment are not 
unexpected. However, studies of control-surface characteristics based 
on potential-flow theory should provide a basis for the study of the 
separated-flow phenomena and provide for a logical grouping of the vari- 
ables to be experimentally investigated. 

The present calculations have been based on two-dimensional aero- 
dynamic coefficients and hence do not take into account the effect of 
aspect latlo. Aspect ratio could have an appreciable effect on this 
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type of oscillation since, for the most part, the oscil l ations are of 
a low frequency and approach the range of stability frequencies. How- 
ever, reference 6 shows that the effect of finite aspect ratio for a 
pitching wing in supersonic flow reduces hut does not eliminate the 
instability, and a reasonable assumption is that a similar effect may 
be found for the subsonic case. 

The influence of additional degrees of freedom was not considered 
in the analysis. Actual configurations will normally have more than 
one degree of freedom, but an analysis based on a slngle-degree-of- 
freedom system may be an easily obtained limit for cases of coupled 
flutter. There may, ho^vrever, be circumstances in which only one degree 
of freedom is of Importance. Cunningham (reference 3) showed in a 
recent paper that a slight relaxing of the condition of infinite stiff- 
ness of another degree of freedom did not appreciably influence the 
flutter speed. 

Another factor which might be of iD5)ortance is the thickness of 
the control surface. The aerodynamic coefficients used here wre based 
on the concept of replacing the -vrlng and control surface by an infinitely 
thin mean camber line. The effect of tMckness requires further 
investigation. 


CONCLUSIONS 


Single -degree-of -freedom-flutter calculations which show the effects 
of various independent parameters, namely, Mach number, location of hinge 
axis, aerodynamic balance, and structural damping are presented for a 
control surface. The following conclusions may be enumerated: 

1 . Calculations based on 'unsteady potential-flow theory indicate 
the existence of slngle-degree-of-freedom flxrtter of a control surface. 

2 . Flutter of a control surface alone is more likely to occur for a 
configuration operating at high siibsonic or low supersonic speeds and at 
high altitudes than at low speeds eind low altitudes. 

3. Structural damping has a beneficial effect, since it raises the 
flutter speed appreciably. The use of struct'ural damping may be a 
convenient method of eliminating single-degree-of-freedom flutter. 

i)-. The unstable oscillation is still possible if the 
is aero dynamic ally balanced. 


control surface 
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5 . Ttiese results haTre several practical limitations . Such, effects 
as separated flow, aspect ratio, coexistence of other degrees of freedom, 
eind finite thickness of control surface were not considered and could he 
significant. However, the results should prove useful as a basis for 
further experimental study. 


Langley Aeronautical Laboratory 

National Advisory Committee fop Aeronautics 
Langley Field, Va., August 17, 1951 
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Figure 1.- Plot of flutter-speed parameter 
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for various axis-of-rotation locations of a control surface 


in incompressible flow. 
















Figure 2.- Continued 
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Figure 3 .- Reduced velocity l/k at whlcli oscillation vlll occtir plotted 
against location of control- surface axis of rotation c for Incom- 
pressible fl,ow. 








Figure 7- - Plot of flutter-epeed parameter — ^ against - 
different amounts of aerodynamic ‘balance. M => Oj c 
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Figure 10 . - Plot of frequency ratio ( against an inertia param- 
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Figure 11.- Asymptotic value of inertia parameter 

nuniber. c =0.7. 


jrpb^ 


against Mach 
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